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I \ Abstract 

■ The proton analyzing power Ay in pion production reaction 'jfd ^He has 

■ been calculated including one- and two-body meson production mechanisms with a 
I proper treatment of the three-nucleon dynamics and an accurate solution of the 3N 
' bound-state problem for phenomenological two-nucleon potentials. In the region 

^ . around the A resonance, the structure of the analyzing power can be understood 

once interference effects among amplitudes describing intermediate AA^ formation in 
different orbital states are considered along with the additional interference with the 
' ^ ■ S-wave pion production amplitudes. Then, the inclusion of three-nucleon dynamics 

O I in the initial state produces the structure of the analyzing power that has been 

^ I observed experimentally. 

! It has been known for almost twenty years that the experimental structure of the 

^ I analyzing power Ay for the 'p'd —>■ tt^ ^B.e reaction (in the region around the A resonance) 
d ' exhibits a rapid variation at the cm. polar angle of Oc.m. = 90°, with the appearence of 
an additional (second) peak and changing the sign as one moves from lower energies to 
the A-resonance region. To our knowledge, this behavior of the proton analyzing power 
has never been interpreted theoretically before, while the understanding of the underlying 
mechanisms is of great importance to describe and understand the basic NN — > NNtt 
inelasticities in the most simple "complex" nuclear environment, the three-nucleon system. 
Single pion production represents the first process that needs to be understood in order 
to consider more complicated hadronic processes that occur at moderately higher energy, 
such as double-pion production, or heavier meson (e.g., rj) production. In addition, single 
pion production could also reveal mechanisms involving multinucleon exchanges, thus 
providing important information that bridges intermediate-energy phenomena to low- 
energy three-nucleon and multinucleon forces, to be used in the study of low-energy nuclear 
systems. 

Theoretically, pion production in proton-deuteron collisions was studied long time ago 
with the deuteron model [1 , which describes the pd — > 7r~^t reaction in terms of the dif- 
ferential cross section for the more elementary pp Tc^d process. The model has been 
subsequently refined by Locher and Weber 0, who introduced additional mechanisms 
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based on the elastic nd nd cross section, and by Fearing ^ who gave phenomenological 
estimates of the effects introduced by optical-model distortions. The deuteron model has 
been later extended by Germond and Wilkin |4j who developed a simple spin-dependent 
analysis in plane waves by which they predicted the deuteron tensor analyzing power 
T20 to approach its geometrical limit — at threshold, based on low-energy theorems. 
Recently, Falk j3] developed further the deuteron model in its spin dependent form and 
performed a phenomenological analysis of polarization observables near threshold. Cal- 
culations including explicitly A dynamics and meson-exchange processes were pioneered 
in the late '70 [Hj, but the complexity involved in the computation of the corresponding 
amplitudes forced the authors to restrict themselves with only exploratory calculations. 
Laget and Lecolley [7j performed calculations that tried to go beyond one- and two-body 
pion production mechanisms and included effective three-body mechanisms via meson 
double rescatterings, showing beneficial effects at higher energies, well beyond the A res- 
onance. Results around the isobar resonance, however, were not quite satisfactory, and 
especially the proton analyzing power Ay, i.e., the asymmetry of the reaction yield for the 
opposite proton spin directions, 

was in poor accord with the data IH|. 

The approach involving the 3A^ dynamics has been pushed forward recently in a series 
of papers (Hj, where explicit A excitation mechanisms were treated in the three-nucleon 
system using the few-body techniques [TO]. Intermediate A propagation and rescattering 
were described phenomenologically using an isobar complex-mass scaling that has been 
tested and parametrized on the pp — >■ 7r"'"c? reaction [llj. Such test calculations have shown 
that, for the description of the proton pionic capture by protons, the interplay between 
(A-mediated) P-wave mechanisms and S-wave two-body rescattering mechanisms (in the 
isoscalar and isovector channels) is a fundamental feature that has to be taken into account 
for the overall reproduction of both spin-averaged and spin-dependent observables. Anti- 
symmetrization prescriptions due to the identity of the nucleons were fully taken into 
account. 

Nowadays, the calculations can be performed with more accurate knowledge of the 
nuclear wave functions (WF), and with including a much larger number of intermediate 
three-nucleon states. Furthermore, the nucleon-deuteron initial-state interaction (ISI) 
can be calculated through, e.g, the Alt-Grassberger-Sandhas (AGS) scheme |T^, with 
realistic nucleon-nucleon transition matrices. The complete calculation of the reaction 
represents a very complicated theoretical task, and involves, due to the variety of the 
few-body intermediate states that have to be expanded in partial waves, computation of 
an exorbitant number of multidimensional integrals. 

The amplitude of the pd —>■ ^He reaction (or the inverse pion absorption process) 
is determined by the matrix element 

Af =^-^ siq^Hm.'Pl) , (2) 

where ipt represents the three-nucleon bound-state and s{^.^''Pd\ refers to a properly 
antisymmetrized interacting nucleon-deuteron state in the initial channel. Diagrammati- 
cally, this amplitude is illustrated by Fig. [T] Within the model discussed here, it includes 
the proton-deuteron ISI, the bound-state two- (deuteron) and three-nucleon (^He) WF, 
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Figure 1: Diagrammatic representation of the pion production amplitude. 



the matrix element of a 3A^ permutation operator that gives rise to the exchange contri- 
bution to the amplitude, and an operator of pion production on a pair of nucleons. Note 
that in formulation of the AGS formalism employed in our calculations, the 3A^ break-up 
process is effectively included in the ISI block in Fig. ^ 

The meson production mechanisms are given by the two-body (Fig. |2t',b,c) and one- 
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Figure 2: Meson production mechanisms considered. 
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body (plain irNN vertex in Fig. |2Ji) contributions. The A excitation mechanism (Fig. |2b) 
has the structure 

-^A = Vnn-naGaF^antt , (3) 

where Vnn-na is the NN — > A^A transition potential generated by vr- and p-exchanges, 
Ga effectively describes the intermediate ANN propagation, and Fan-k is the AA^vr tran- 
sition vertex function. All the corresponding details including the extended formulas 
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representing the mechanisms in the three-nucleon partial waves have been pubhshed else- 
where jni CD] • The asymptotic pion wave is not decomposed in partial waves but treated 
directly in three dimensions. An important quantity for the present calculation is the 
maximum orbital momentum /a in the AA^ subsystem. We found that /a = is not 
sufficient to describe polarization observables, and one has to go up to /a = 2, at least. 

In the A resonance region, the spin structure of the reaction amplitude is depen- 
dent strongly also on the interference effects with the other mechanisms, particularly 
the S-wave two-body pion production contributions shown in Fig. ^i,h. Each of these 
mechanisms can be schematically described as 

Att = Fnn-kGt,\n-k-Ntt + time reordering , (4) 

where Fattvtt is the vrA^A^ vertex, G-,, is the pion propagator, and AAr-Tr-ArTr is the correspond- 
ing contribution to the half-off- shell ttN scattering amplitude. By "time reordering" we 
denote a term coming from "backward in time" propagation of the intermediate-state 
pion , in which case Aat^-attt is replaced by a — 2t:N vertex amplitude for double-pion 
production. 

The starting point for the construction of AAfTr-AfTr is the phenomenological low-energy 
(5- wave) pion-nucleon interaction Hamiltonian ^3] 

Hint = ^i'l^Ti^ ■ W X d^ct>] + [0 ■ 0] . (5) 

The calculations herein illustrated have been performed with the set of parameters 
used in the previous study [H]. Namely, for the isovector coupling constant in Eq. (jSj), 
we have Av'/47r = 0.045, whereas the vrA^A^ and vrA^A coupling constants are equal to 
f^^j^/An = 0.0735 and /^ata/^tt = 0.32, respectively. Our present treatment is based on 
Ref. ^3]. The isospin-odd contribution is given in terms of a p-exchange model, while 
the isospin-even term is described as the combined effect of a phenomenological short- 
range (SR) process (Fig. ^) and an effective scalar-meson (a) exchange of Fig. |2t. The 
two latter mechanisms act in opposite directions and almost cancel each other in case 
of the free pion-nucleon scattering. However, an off-shell enhancement of the probability 
amplitude takes place in the scalar-isoscalar channel for pion production off a bound 
system. A question still remains open, if this form mimics other physical effects, such 
as the 2N SR-exchange contributions proposed in Ref. [1^] (see also Ref. [16j, where the 
current status of theory and phenomenology of pion production in A^A^ collisions close to 
threshold is reviewed). 

Further details (in particular, the treatment of the fully antisymmetrized matrix ele- 
ments with respect to the Faddeev three-nucleon wave function) are not discussed here 
and can be found in the papers [3 El E]- The three-body dynamics in the initial 
state has been described following the treatment [T7j of the low-energy proton-deuteron 
scattering, with the set of the AGS equations being reduced through application of the 
spline-interpolation technique to a linear equation system. 

In Fig. 01 the results for the proton analyzing power are compared with what has 
been observed in experiments. For the three different phenomenological nucleon-nucleon 
potentials, one observes an additional peak at Oc.m. — 100° when going from i?iab = 350 
MeV towards higher energies. This behaviour is fully explained by the model calculation, 
once the ISI involving the three-nucleon dynamics is taken into account. Otherwise, if 
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Figure 3: Proton analyzing power for three different proton lab energies. Solid, dash- 
dotted and dotted lines denote our complete calculations for Bonn-B, Bonn- A, and Paris 
A^A^ potentials, respectively, while thin solid line corresponds to the plane-wave calculation 
(without the pd 1ST). Full triangles, blank circles, blank triangles and full circles mark data 
of Refs. [HllIHllini and [201, correspondingly. 
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one neglects these three-body effects in the initial channel, the structure of the observable 
is poorly reproduced. 

In generation of this structure, another quantity that plays an important role is the 
angular momentum of the intermediate AN state. While /a = is sufficient in description 
of the general behaviour of the excitation function (see Ref. jHI), one must consider higher 
AA^ partial waves, at least up to /a = 2, in case of the analyzing power Ay. As follows 
from Fig. |31 there is no need, however, to go beyond this value, since the numerical 
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Figure 4: Results for different values of the maximum AA^ orbital momentum. Dotted, 
solid and dashed-dotted curve refer to calculations with /a = 0, 2, and 4, respectively, 
performed with the Bonn-B potential and inclusion of the 3A^ ISI. 

convergence is already reached with the AA^ D-waves. 

Apart from the analyzing power Ay, also the angular differential cross-section has 
always been a challenge for the theory. Recent COSY data [21] provided new informa- 
tion about the shape of the cross- sect ion, being at the same time consistent with earlier 
TRIUMF measurements jH]. AH these data compare fairly well with our calculations, as 
shown in Fig. |S| Indeed, almost all of the data points fall in the strip determined by the 
calculations performed with different A^A^ potentials. It is worth noting that for higher 
energies and backward angles, the data seem to favor more the Paris potential, while in 
the lower part of the energy range considered, the Bonn potentials (Bonn-B, in particular) 
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Figure 5: Differential cross-section of the pd —>■ 7r° ^He reaction. Data of Refs. [21] and [S] 
are shown by full circles and triangles, respectively. Notations for the curves are the same 
as in Fig. |21 
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give rise to the results which are closer to the data. 

To summarize, we have described the pionic capture processes in nucleon-deuteron 
collisions employing elementary production mechanisms obtained from the non-relativistic 
reduction of the phenomenological ttN and ttA Hamiltonians. The corresponding matrix 
elements between the three-nucleon states have been evaluated within a large-basis space, 
consisting of 464 three-nucleon partial waves. Total angular momenta of the system up 
to 7/2 for both the parities have been considered. The prescriptions required by the 
Pauli principle have been taken into account through the application of the permutation 
formalism to the three-nucleon system. Finally, these matrix elements have been folded 
with the nuclear WF to obtain the probability amplitudes, with the complete details of 
the procedure being outlined in Refs. CHI HI]- In the calculation of the plane- wave 
part of the amplitude, we mainly follow these papers. The only difference is modification 
of the standard Jacobi momentum set for the system through the replacement of the 
pion mass by its total energy, what is more adequate to the kinematics of the present 
calculations. An extended paper explaining this and some other ingredients of the model 
is in progress. However, much more important turns out to be our regard of the 3A^ 
dynamics in the initial proton-deuteron state, which was not extensively investigated in 
the previous treatment of the reaction. 

The results obtained suggest that one achieves a fairly good reproduction of the phe- 
nomenology, when including, in addition to the A-rescattering mechanism, the p-exchange 
process in the S-wave, as well as another S'-wave mechanism that generates a significant 
amount of interaction in the scalar-isoscalar vrA^ channel. On the contrary, the one-body 
term (the off-energy- shell ttNN vertex) adds contributions that are significantly smaller. 
For the scalar-isoscalar component we have considered an off-shell model ^3], which has 
been used before. We do not claim here that this particular way to generate the flux 
in the scalar-isoscalar channel is more realistic than other alternatives (e.g., the above- 
mentioned SR-exchange contributions proposed in Ref. A study on how really the 
flux is generated in the scalar-isoscalar channel is beyond the scope of the present work, 
but can represent at this point a subject for future investigations. We hope that studies of 
spin observables close to threshold of this specific reaction can shed light on that problem. 

This work shows that, with a combination of two-body mechanisms, one can explain 
the structure of the proton analyzing power Ay, which has been accurately measured 
about 20 years ago and represented since then a challenge for theorists, being a difficult 
observable to interpret. We conclude that it is possible to reproduce qualitatively the 
nontrivial energy dependence of Ay around the A resonance, provided that the three- 
body nuclear dynamics in the initial state is taken into account. It is worth mentioning 
also that we have made use of the techniques that have been developed, in fact, mainly for 
low-energy few-nucleon applications. It is remarkable that they can be successfully used in 
this energy regime, where the use of the three-body non-relativistic quantum-mechanical 
equations and the concept of the standard A^A^ potentials is, in some sense, questionable. 

L.G.L. is indebted to the University of Padova and INFN for their kind hospitality 
and support in October-November, 2003 and April-August, 2004. 
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